Effect of the In-Plane Magnetic Field on Conduction of the Si-inversion Layer: 

Magnetic Field Driven Disorder 
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We compare the effects of temperature, disorder and parallel magnetic field on the strength of the 
metallic-like temperature dependence of the resistivity. We found a similarity between the effects of 
disorder and parallel field: the parallel field weakens the metallic-like conduction in high mobility 
samples and make it similar to that for low-mobility samples. We found a smooth continuous effect 
of the in-plane field on conduction, without any threshold. While conduction remains non-activated, 
the parallel magnetic field restores the same resistivity value as the high temperature does. This 
matching sets substantial constraints on the choice of the theoretical models developed to explain 
the mechanism of the metallic conduction and parallel field magnetoresistance in 2D carrier systems. 
We demonstrate that the data for magneto- and temperature dependence of the resistivity of Si- 
MOS samples in parallel field may be well described by a simple mechanism of the magnetic field 
dependent disorder. 

PACS: 71.30.+h, 73.40.Hm, 73.40.Qv 



I. INTRODUCTION 

The apparent metallic-like temperature dependence of 
the resistivity is observed in various two-dimensional 
(2D) carrier systems and remains a focus of a broad 
interest because it challenges the conventional theory of 
the metallic conduction. The effect manifests itself in 
a strong drop of the resistivity over a limited range of 
temperatures, T — (0.5 — 0.05)Tp, from a high temper- 
ature value (call it 'Drude' value) ptugh = p{T ~ Tp) to 
a low temperature one po = p{T ;$ 0.05Tp), here Tp is 
the Fermi temperature. Upon lowering the temperature 
further, a strong 'metallic' drop in p was found to cross 
over to the conventional weak localization type depen- 
dence p|-|T3|. Recently, it was demonstrated p| , |l3| that 
the 'metallic' drop is not related to quantum interference 
and in this sense should have a semiclassical origin. 

Magnetic field applied in the plane of the 2D system 
causes a dramatic increase in resistivity [jTjl-jbSl ■ It was 
proven experimentally jl^,^o| that the magnetoresistance 
(MR) in Si-MOS samples is mainly caused by spin-effects, 
though certain contribution of orbital effects pl| is no- 
ticeable at very large fields [^(|, bigger than the field of 
the spin polarization. The two effects, strong metallic 
drop in resistivity and parallel field magnetoresistance 
remain puzzling. 

In this paper we demonstrate that (i) the strength of 
both effects is well matched and (ii) the action of mag- 
netic field on conduction is similar to that of disorder and 
to some extent to that of temperature. In particular, in 
the presence of in-plane magnetic field the strong tem- 
perature dependence of the resistivity for a high mobility 
sample transforms into a weaker one and shifts to higher 
temperatures; both these features are typical for low mo- 
bility samples. Next, (iii) the effect of the magnetic field 
on the resistivity is continuous and shows no signatures 
of a threshold. Finally, (iv) increase in either temper- 



ature or parallel magnetic field restores the same high- 
temperature ('Drude') resistivity value while transport 
remains non-activated. Matching of the actions of the 
above controlling parameters sets substantial constraints 
on the choice of theoretical models. 

We performed measurements on a number of different 
mobility (100) n— Si-MOS samples, because this system 
demonstrates the most dramatic appearance of the dis- 
cussed effects. Table | below shows the relevant parame- 
ters for the four most intensively studied samples whose 
mobilities differ by a factor of 25; parameters for other 
samples were reported in Refs. [|||,|l^,^|. 

TABLE I. Relevant parameters of the five samples. 
Ppeak [m 2 /Vs] is the peak mobility at T = 0.3 K, p c is 
in units of h/e 2 , and n c is in 10 11 cm -2 . 



sample 


/^pcak 


n c 


Pc 


Si 9 


4.3 


0.75 


3.2 


Si 15 


4.0 


0.88 


2 


Si 43 


1.96 


1.4 


0.67 


Si 39 


0.51 


3.5 


0.32 


Si 46 


0.15 


8.3 


0.3 



In Figures 1 a, c we plotted the temperature depen- 
dence of resistivity p for the highest and lowest mobil- 
ity samples. At high carrier density (the lower part of 
Fig. 1 a) the resistivity decreases and demonstrates the 
metallic-like behaviour. The magnitude of the drop in 
resistivity, (/Chigh — Po)/ POt obviously depends on the dis- 
order and varies from a factor of 5 (the lowest curves in 
Fig. la) to a few % in Fig. 1 c. As density decreases 
below a "critical value" n c , the character of the tem- 
perature dependence changes to the insulating one. The 
changes in the character of p(T) at n = n c are reminis- 
cent of a typical metal-insulator transition (MIT) in 3D 
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systems which is not expected to occur in 2D system. In 
agreement with earlier observations [^p|JT^| the increase 
in disorder causes the drop in resistivity to weaken, the 
'critical' density n c to increase and the 'critical' resistiv- 
ity p c to decrease. 
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FIG. 1. T— dependence of the resistivity (a) for high 
mobility sample Sil5 at B = 0, and (b) at Bu — 12 T; 
(c) for low-mobility sample Si46 at B = 0. Inset shows 
p(n) for T — 1, 2, 4, 8, 16, 32, 62 K (arrows indicate direc- 
tion of the growth in temperature). Horizontal dashed 
lines correspond to a 'critical density' n c . Dash-dotted 
line in the panel (a) marks T = O.lTp. The densities (in 
units of 10 11 cm -2 ) on the panel (a) are 0,76, 0.79, 0.85, 
0.89, 0.90, 0.91, 0.96, 1.01, 1.07, 1.12, 1.18, 1.23, 1.29, 
1.40, 1.51, 1.62, 1.73, 1.84, 1.95; (b) 0.95, 1.07, 1.18, 1.29, 
1.51, 1.73, 2.06, 2.28, 2.61, 2.94, 3.39, 3.72; and on the 
panel (c) 3.85, 4.13, 4.83, 5.53, 6.23, 7.63, 9.03, 10.4, 11.8, 
13.2, 16.0, 18.8, 21.6, 24.4, 30.0, 37.0. 

However, in contrast to the common believe that only 
high mobility samples are subject to the MIT, we demon- 
strate in the inset to Fig. 1 c that even for such low 
mobility sample as Si46, there is a clear crossing point, 
n = n c , where dp(T)/dT changes sign. The difference 
from that for high mobility samples (Fig. 1 a) is the shift 
of the events to much higher temperatures, to the range 
T = 1 — 60K , whereas for lower temperatures T < 1 K 



the weak localization upturn in the resistivity sets in. 
This upturn is a subject of more detailed consideration 
p3[ and is not discussed here. The increase in the tem- 
perature scale by about 10 times as a function of disorder 
is anticipated since the Fermi energy, Ep, at the critical 
carrier density for disordered sample Si46 are 10 times 
bigger than that for Sil5; correspondingly, in both sam- 
ples the MIT is observed in similar ranges of (T/Ep). 
Comparing Figs, la and lc we arrive at the conclusion 
that the overall temperature dependence of the resistiv- 
ity for Si-MOS samples may be roughly scaled by taking 
into account disorder dependence of the critical density 
n c and the corresponding energy scale. We ignore at the 
moment a disorder dependence of the strength in the re- 
sistivity drop, which will be considered further. 

II. THE EFFECT OF THE PARALLEL 
MAGNETIC FIELD ON CONDUCTION 

A. Magnetic field driven Mi-transition? 

It was found in Refs. Jl7],[l^,^| that for a fixed car- 
rier density an increase in Bn causes dp/dT to change 
sign. At first sight, this effect is reminiscent of a magnetic 
field driven transition; an analogy to the superconductor- 
insulator transition was discussed in Ref. and the cor- 
responding field value is often called "critical magnetic 
field" B c . In Figure 2 we reproduce similar behavior for 
n— Si-MOS samples. 
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FIG. 2. a, b) Resistivity traces vs Bu for different tem- 
peratures and for two densities. Sample Si43. Vertical ar- 
rows mark interval of fields which may be assigned to the 
'critical' field value B c . c) B c vs carrier density. Shadowed 
corridor corresponds to the interval of crossing. Dashed 
line is a linear fit. 
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Similar to that reported in Refs. @J||§7|], the p{B) 
traces for different temperatures cross each other around 
a 'critical' magnetic field, B c , though in an extended in- 
terval rather than in a single point. As carrier density 
ri and magnetic field increase, the interval of crossing 
broadens and finally becomes infinite for B > 6 T. Den- 
sity dependence of B c is plotted in Fig. 2 c in the range 
of fields where the crossing interval is finite. 

In order to understand the driving mechanism of the 
magnetoresistance we performed more detailed studies 
of the magnetoresistance and measured its temperature 
dependence at various values of the in-plane field. 



'critical' density, n c , separates regimes of temperature 
activated (i.e. 'insulating', n < n c ) and non-activated 
('metallic', n > n c ) conduction p<| . 

The figures clearly demonstrate again that the paral- 
lel field progressively increases the critical carrier den- 
sity n c and decreases the corresponding critical resistiv- 
ity value p c . For magnetic fields Bu > 6T (not shown), 
the crossing point spreads into a density interval, though 
still shifting to higher densities. This is typical for low 
mobility samples and typical for the effect of the disorder. 



B. Analogy between the effects of parallel field and 
disorder on conduction 

It is rather instructive to trace in detail the influence of 
the magnetic field on the density and temperature depen- 
dences of the resistivity. The most striking and known 
result of the in-plane magnetic field is the increase in 
the sample resistivity. Figure 3 demonstrates that the 
mobility degrades in magnetic field smoothly, without a 
threshold. Besides decreasing sample mobility, the mag- 
netic field causes its maximum to shift to higher density. 
Both these effects are typical for the B = case when 
the mobility is varied by changing the disorder; at low 
carrier density this is usually described in terms of the 
increase in the number of scatterers, rij oc 1/ fjP cak p8] . 
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FIG. 3. Mobility vs carrier density for a fixed temper- 
ature T = 0.29 K and for different magnetic fields B\\ (in- 
dicated in the figure). 

In Figures 4, we show the resistivity traces vs density 
for four fixed temperatures in the range 0.3 to IK. All 
four groups of curves for different B\\ values demonstrate 
a critical behaviour. The well defined crossing point or 
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FIG. 4. Resistivity vs carrier density for 4 fixed 
temperatures (indicated on the figure) measured at 
B\\ = 0, 1, 2 and 3 Tesla. Arrows at each panel mark 
the critical density and resistivity. 



From Figure 4, we determined the coordinates of the 
crossing point, n c and p c , for each Bu value and plotted 
them in Figs. 5. These two curves n c (B) and p c (B) have 
a transparent meaning: they separate the regimes of the 
temperature activated and 'non-activated' conduction on 
the planes p(B, n) We find again that both field de- 
pendences are qualitatively similar to the ones known for 
the action of disorder (e.g., inverse sample peak mobility 
l/fi, or 1/M value) 
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FIG. 5. Critical resistivity (a) and critical carrier den- 
sity (b) vs parallel magnetic field. Full line is a linear fit 
with adjustable slope and with n c = 0.883 x 10 11 cm' 2 . 
Dashed line is a guide to the eye. 

Figure 6 shows for comparison the p c (n c ) dependence 
measured with the sample Sil5 in various parallel fields 
(calculated from Fig. 4) and p c {n c ) reproduced from 
Ref. H for different samples at B = 0. A clear anal- 
ogy is seen between the decrease in p c caused by disorder 
and magnetic field. 
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FIG. 6. "Critical" resistivity value vs "critical" den- 
sity. Closed squares show the dependence driven by mag- 
netic field, open squares are for the one driven by disorder 
(plotted for 9 different samples, based on the data from 
Ref. |^]). Line is the guide to the eye. 

Finally, in order to illustrate the conclusion that the 
action of the magnetic field is similar to that of the dis- 
order, we plotted in Fig. 1 b, for comparison, p(T) for 
the same high mobility sample Sil5 but measured with 
applied parallel field of 12 T. In light of the above discus- 



sion, the transformation of Fig. 1 a into Fig. 1 b can thus 
be treated as a result of a progressive increase in electron 
scattering introduced by the parallel magnetic field, in 
steps like those shown in Figs. 4 and 3. The similarity 
between Figs. 1 b and 1 c supports this conclusion. 



C. Magnetoresistance in parallel field and the 
o-factor 



A number of experimental studies (e.g., Refs. |19| , [20[ ) 
demonstrate that the parallel field magnetoresistance in 
Si-MOS samples is a spin rather than orbital effect. 
There were, consequently, a number of empirical at- 
tempts to extract the g*-factor value from different fea- 
tures of the magnetoresistance R(B^). For comparison, 
we present in Fig. |?] the density dependences of two char- 
acteristic magnetic fields: B sat , corresponding to the sat- 
uration of R(B\\) (for Sil5 and Si43) from Ref. @, and 
the 'critical magnetic field' B c (for Sil5) re-plotted from 
Fig. 5. 

First of all, it is clear, that B sat does not depen- 
dent solely on the carrier density, but does also depend 
on the disorder. Whereas the empirical dependences 
B sa t — s(n — rib) for different samples have almost the 
same slope s, the offset rib was found to be inversely 
proportional to the sample mobility ppj ]. For different 
samples, rib varies from 0.5n c for the highest mobility 
sample Si9 to 0.7n c for Sil5, 0.8n c for Sil2, and 0.96n c 
for Si43. On the other hand, the 'critical field' B c , by 
definition, extrapolates to zero at n — n c (B c ). For this 
reason, the two curves, B sa t(n) and B c (n) usually inter- 
sect. Obviously, rib cannot be identified with n c . At low 
densities, -B sa t exceeds B c and the saturation of the mag- 
netoresistance takes place in the insulating regime. At 
higher densities, the saturation occurs in the 'metallic' 
regime; the borderline for the sample Sil5 corresponds 
to n = 1.75 x 10 n cm- 2 . 

For comparison, the dashed-dotted line in Fig. fj] repre- 
sents the calculated density dependence of the field corre- 
sponding to the complete spin polarization of conducting 
electrons: 



B„ i = 



2E F 

9*vb 



g*rrr 



2 



(1) 



where the renormalized <?*m*-values for conducting elec- 
trons were directly measured in Ref. |50| as a function of 
the carrier density. 

For a given density n, B sat {n) is usually less than the 
field of complete spin-polarization B po \ of the conducting 
electrons. The deficit, B po \ — B satl increases as sample 
mobility decreases for example, B po \ — £> sat w 4T 
for Si43. Whereas i? po i is solely determined by the mo- 
bile carrier density, B sat turns out to be also dependent 
on disorder. 



4 



For high mobility samples, over a certain range of den- 
sities (which is (1.2 - 2.2) x 10 n cm" 2 for Si 15) the two 
quantities, -B po i and B sat , are rather close to each other. 
This coincidence is in a good agreement with the obser- 
vation by Vitkalov et al. pl| who found that in the range 
of densities n — (1.54 — 4.5) x 10 11 cm~ 2 , the frequency of 
the Shubnikov-de Haas oscillations in tilted fields doubles 
at a field which is equal to B sat to within 5%. However, 
for more disordered samples (such as e.g. Si43 in Fig. |^), 
the saturation occurs in magnetic fields essentially lower 
than the polarization field i? po i, because of the twice as 
large n& value (data on rib for a number of samples may 
be found in Ref. We believe therefore that the coin- 

cidence of i? sa t and B po \ for some high mobility samples 
and in a limited density range is rather occasional. 

In the analysis of the £? sat (n)-data in Ref. pi]] , it was 
assumed that B aat remains equal to B po i in the n — > n c 
limit. With this assumption authors arrived at the con- 
clusion that g* diverges as the density decreases, and 
that a ferromagnetic transition takes place at the MIT. 
As Fig. 7 shows, the coincidence of B sat and B po i fails for 
densities lower than 1.3 x 10 n cm~ 2 (even for the high 
mobility sample). The <?*-factor for conducting electrons 
measured in Ref. (3TJ| from Shubnikov-de Haas effect does 
not diverge at n c , but gradually grows as the density de- 
creases; the growth is anticipated within the Fermi-liquid 
theory p5|,p(|. 




n (10 cm ) 

FIG. 7. B c and i3 sa t for the sample Sil5 (full diamonds) 
and -B sat for Si43 (empty squares) measured in the range 
B\\ = 1.3 — 13 T and extrapolated schematically over a 
wider range of B\\. B po i is the field for complete spin po- 
larization of the conducting electrons, calculated according 
to Eq. (Jl|) and using g*m*{n) values measured in Ref. | pp| . 

One might conclude that the non-zero value of rib 
is the only reason for the difference between -B sat and 
_Bpoi, and that by subtracting rife the two parameters 
may be brought into agreement. However, such proce- 
dure performed in Ref. pfl for different samples, has 



revealed that the two parameters are indeed related to 
different physics: the (<7to)mr values estimated from the 
magnetoresistance either as (/i/e)n/_B sat , or (h/e){n — 
nb)/B aa x disagree with the (g*m*)sdH-values measured 
from Shubnikov-de Haas effect. This is also seen in Fig. 
where (i) at low densities the difference between the mea- 
sured B sat (n) values and the B po i{n) curve can not be 
eliminated by a horizontal shift (i.e., by varying rib) and 
(ii) the slopes of the B sat (n) and B po i (n) curves are dif- 
ferent. 

We now consider another characteristic field, B c . 
In Ref. p7d , the £? c (n c )-dependence for 2D holes in 
(311)GaAs was noted to correspond to such spin align- 
ment in the 2D carrier system where the spin-minority 
population drops below a threshold (which is approxi- 
mately independent of the total 2D hole density or the 
magnetic field); the latter is of the order of the critical 
density at zero field n c (B = 0) ^f\ . In order to verify 
such possibility for Si-MOS samples, we fitted in Fig. 5 b 
the n c (B) data for the high mobility sample Sil5 with a 
linear dependence n c (B) — rib = cB, with two adjustable 
parameters rib and c. We find that the offset, rib, is in- 
deed equal to the anticipated value n c {B = 0) = n c (for 
sample parameters, see Table 1), however, the slope c = 
0.13 x 10 n cm _2 /T, disagrees with such interpretation. 
With this interpretation, the slope c = dn c /dB c would 
be equal to (g*m*)g v e/2h (where g v — 2 for n-(lOO)-Si); 
we obtain then (1/2)(<?*to*)mr = 0.27 over the displayed 
range of densities n = (0.9 — 1.5) x 10 n cm~ 2 . We re- 
peated the same procedure with the data from Fig. 2 b 
and found, correspondingly, (l/2)(g*m*)MR ~ 0.48 for 
the range n = (2.2 - 2.7) x 10 11 cm -2 . 

The increase in the slope, dB c /dn with carrier den- 
sity would mean the decrease of g*m* for lower densities. 
Such decrease contradicts the results of direct measure- 
ments of the (g*m*)sdn for conducting electrons from 
Shubnikov-de Haas effect |l7|j3C|l ; the contradiction in- 
dicates that the slope of the B c (n c )— dependence in Si- 
MOS samples might have a more complex interpretation 

Earlier ©] we concluded that the saturation of the 
magnetoresistance (MR) for Si-MOS samples in paral- 
lel field _B sat is more related to the g*-factor of localized 
electrons than that of the conducting ones. We think 
the same is true for the magnetoresistance at B = B c . 
To illustrate this conclusion, we compare in Table II the 
g*m* values (i) directly measured for mobile electrons 
from Shubnikov-de Haas effect (ii) the data which 
follow from the slope, (dn c /dB c ), of the curves in Figs. 2 
and 5, and (iii) the ones which follow from the slope 
(dn/ dB sa t) of the density dependence of the saturation 
field g(]]. 

At r s = 5.3, the g*m* values determined from 
Shubnikov-de Haas oscillations and from MR in paral- 
lel field are comparable. However, at r s = 7.6 they differ 
by a factor of two: g*m* determined from the SdH effect 
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increases with r s as expected for conducting electrons 
p0| ]; in contrast, (<7*m*)B c and (<7*w*)B sat are close to 
the value 2m\, (mj = 0.19 is the bare band mass) antici- 
pated for the localized electrons. 



TABLE II. Comparison of the (g*m*) and g* measured 
from Shubnikov-de Haas effect (labeled 'SdH') |3(J with 
values estimated from the MR in parallel field. Label l B c ' 
denotes the data derived from dB c I dn c , and label '-Bsat' 
denotes the data from dB Ba , t /dn [|20|. 



r s 


i(g*m*)sdH 


i(g*m*) Bc 


i(g*m*) flsat 


5.3 
7.6 


0.44 
0.665 


0.48 
0.27 


0.35 
0.34 



We close this section with a note that the attempts 
of describing B sat and B c in terms of the diverging g*- 
factor |33],[5l) , face with problems of explaining (a) why 
the 'divergence' takes place at progressively (and essen- 
tially) lower r s values (i.e. higher n c values) as disorder 
increases and (b) why the diverging g* values are different 
from those measured in Shubnikov-de Haas effect. The 
comparison presented in this section illustrates a rather 
complex behaviour of the magnetoresistance in the paral- 
lel field; definitely a thorough theoretical consideration 
is needed in order to use this effect for extracting the 
electron spin properties. 



D. Comparing the effects of and temperature on 
conduction 



We would like to point to the interesting similarity be- 
tween the action of the parallel magnetic field and the 
temperature on the resistivity. Figure || demonstrates 
that for sufficiently high carrier density, n 3> n c , the ap- 
plication of high magnetic field increases resistivity by 
about 5-6 times and restores nearly the same zero-field 
'Drude' resistivity value as the temperature does: i.e. 
p(B > B sat ,T < T F ) = p(B = 0,T > T F ). The sim- 
ilarity between the effect of the field and temperature 
holds only while conduction remains 'non-activated' (cf. 
Figs. 5). 

The dashed line in Fig. ||b, depicts the boundary (same 
as in Fig. 4) between the temperature-activated and 'non- 
activated' transport regimes. When the carrier density 
becomes lower than n c (i? max ) but is still bigger than 
n c (B — 0), the parallel field drives resistivity through 
this boundary at some field B c (n). As a result, the con- 
duction mechanism in magnetic field changes and the 
similarity between Figs. ||a and ||b breaks. For this rea- 
son, the uppermost curves in Fig. ^b show much bigger 
increase in the resistance with field than with temper- 
ature 1 32 1 . It is noteworthy, that at low density, the 



regime, at B sa ± > B c , whereas at higher density the sat- 
uration occurs in the 'non-activated' metallic regime (see 
Fig. @. 
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FIG. 8. Resistivity vs temperature (left panel) and 
parallel field (right panel) for different densities. Dashed 
line indicates the boundary for activated/nonactivated 
regime. On the left panel the densities are, from bot- 
tom to top: 1.84, 1.62, 1.40, 1.29, 1.23, 1.18, 1.12, 1.07, 
1.01, 0.96, 0.91, 0.90, 0.89, 0.85, 0.79, 0.76 x lO^cm" 2 . 

For the same reason, a very low magnetic field may be 
sufficient to drive the conduction from 'non-activated' to 
the 'activated' regime if the density n is chosen slightly 
less than n c . In other words, this may be observed when 
p(B = 0) is chosen very close to but less than p c (B = 0) 
(see Figs. 5 a and 8 b). In light of this note the destruc- 
tion of the metallic conduction in 'arbitrary small field' 
|p"5f and the (Bn/T) scaling of the resistivity at high p 
reported in Ref. ||l5f obtain a transparent interpretation 



E. Magnetic field driven disorder 



magnetoresistance saturation takes place in the activated 



It was verified experimentally |1S|J20[| that in Si-MOS 
structures the parallel field couples to the electron spins 
only. Our results, therefore relate the action of the spin 
polarization and the increase in disorder. This may occur 
for example (i) via weakening of the screening in mag- 
netic field p5| ], (ii) due to inter spin-subband scattering 
p6|- p|p7| , or (iii) via increase in the Coulomb scattering 
caused by 'undressing' (i.e. depopulation and charging) 
of filled interface traps in any model of the interface traps 

As we demonstrate further, with the only assumption 
on the disordering role of the parallel magnetic field, we 
obtain a qualitative description of all presented data. 
We use this assumption as a phenomenological model; 
although we are unable to distinguish which one (or a 
combination) of the microscopic mechanisms is responsi- 
ble for the parallel field magnetoresistance, we evaluate 
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various options in comparison with the data. 

The similarity between the p(T)- and p(-B||)-traces in 
Fig. H seems easy to explain with the mechanism of 
screening which is temperature- p3|— and magnetic 
field-dependent p5| ]. However, comparing this mecha- 
nism with the data we find a number of principle dis- 
agreements: 

• a) The magnitude of the magnetoresistance (factor 
of 5 - 5.5 in Fig. 7 in the 'metallic' range) is larger 
than the screening mechanism can provide ]35|] . 

• b) At sufficiently high parallel field, B > -B Ba t, 
magnetoresistance almost saturates; this field was 
shown to be very close to the complete spin polar- 
ization of 2D carriers |l{J . For the screening mecha- 
nism, even though the spin system is polarized, the 
temperature dependence of resistivity should be as 
strong as in zero field. The dashed line in Fig. 1 b 
corresponds to the density n(B sa t) at which the 
magnetoresistance for this sample saturates. Obvi- 
ously, p(T) depends on temperature much weaker 
than that in Fig. 1 a for B = 0. 

• c) The Hall resistance data [Q also does not sup- 
port the screening mechanism. 

Last but not least: in Ref. p^] a deviation of the weak 
field Hall resistance from its classical value, was observed 
close to the critical density, particularly, on the 'metal- 
lic' side. This means that for n close to n c there arises 
an excess or deficit of delocalized carriers; this observa- 
tion obviously supports the latter scenario (iii). For this 
mechanism to work, the filled traps (or localized carri- 
ers, which is the same for us) should lift in energy with 
parallel field. This may take place for example, if the 
relevant traps form the upper Hubbard band Ml], or if 
the traps are easy or spontaneously spin-polarized. We 
considered the latter possibility in Ref. BOj in order to 
explain the disorder dependence of the field for saturation 
of the magnetoresistance. We note that the spin polariza- 
tion of localized states is favored by the broken inversion 
symmetry of the interface. 

Interface defect charges originating from the lack of 
stoichiometry are intrinsic to Si/Si02 system; their typ- 
ical density is 10 12 cm~ 2 for a thermally grown dioxide 
p6| . As parallel field increases, the band (or upper band) 
of localized carriers should lift in energy and gets 'un- 
dressed' after passing through the Fermi energy. This 
will cause 'turning-on' the charged scatterers and cor- 
responding increase in the scattering rate of the mobile 
carriers. 

The carriers released from the traps will join the con- 
duction band and may be detected via a deviation in the 
Hall voltage as n approaches n c (B); the latter may be 
driven either by decrease in n or by increase in Bu . The 
trapped charge apriory may be of arbitrary sign which de- 
pends on the interface chemistry and growing processing. 



Consequently, the deviation in the Hall voltage may be 
of arbitrary sign as well. Whereas the contributions into 
the Hall voltage of the released electrons and holes may 
compensate each other, it is not the case for the spatially 
separated charged scatterers. The number of 'turned- 
on' scatterers N is therefore expected to be much bigger 
than the deviation in the Hall resistance. In Ref. [ p2[ , 
we observed dVulVn of the order of (1 — 10)% in differ- 
ent samples for the carrier density n ~ 10 11 cm -2 ; from 
this figure a lower estimate for the number of 'turned-on' 
scatterers, Ni > (10 9 - 10 10 ) cm~ 2 , follows. 

It is noteworthy, the inter (spin)-subband scattering 
p6| |s|,|7) is not an alternative option but may be supple- 
mentary to the mechanism of the magnetic field driven 
disorder and charged traps. However, recent experiments 
p0| , [5r| have shown that the scattering times in different 
spin and valley-subbands in (lOO)-Si are close to each 
other even though the degree of polarization is ~ 30%. 
This sets substantial constraints on those two-band mod- 
els which imply the intra-subband scattering times to be 
different. 

For completeness, we mention that the parallel field 
magnetoresistance in Si-MOS samples was measured ear- 
lier by Bishop, Dynes and Tsui 47 and by Burdis and 
Dean [Q (though as a weaker effect on more disordered 
samples) and associated with electron-electron interac- 
tion and Zeeman splitting ]40j] . In this interpretation, 
however, the values of the interaction constant F were 
found to be unomalously large and to decrease as tem- 
perature decreases [H which is not consistent with the- 
oretical expectations. We obtained qualitatively similar 
results on high mobility samples |5^,|5l|. In other words, 
p(i?ll , T) in the 'metallic' range does not scale as a func- 
tion of (B/T) which is expected for the Zeeman term 
5ct = -~0.Q84(F/2)(g*p B B/k B T) 2 @. It seems there- 
fore unlikely the strong magnetoresistance in Si-MOS 
samples to be caused by electron-electron interaction. 

We limited our consideration by strong changes in re- 
sistivity, 5p(B)/p ~ 1. Parallel magnetic field causes 
an additional small effect on the quantum corrections to 
the conductivity, by decreasing the phase breaking length 
]5l| ; the latter though does not contradict our concept of 
the magnetic field driven disorder. 

III. IMPLEMENTATION OF THE MODEL FOR 
ANALYSIS OF OTHER DATA 



In order to verify the above model we compared it with 
all available data for Si-MOS sample and found rather 
good agreement. For the sake of shortness, we present 
here the comparison only with those data for which it 
was explicitly stated in Ref. |34|] that "the enormous re- 
sponse observed at low temperatures is a consequence of 
effects other than parallel magnetic field-induced changes 
in carrier density or disorder strength" . 
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A. Example 1. 

We begin with results obtained by Shashkin et al. 
and reproduce them from Fig. 9 of Ref. p4j. These 
data shown in Fig. |9| "dramatically demonstrate an ex- 
treme sensitivity of the resistance to parallel field" at 
T = 30 mK. The localized behaviour which appears to 
be absent for the curve at B = is restored in a mag- 
netic field" (cited from Ref. |p^| ). Based on this contrast 
it was stated in Ref. H] that the effect of the parallel 
field is something different from the increase in disorder. 

We plot in Fig. || the model curves using a conventional 
temperature activated dependence 



p(T) = p + Piexp(A/T). 



(2) 



The effect of the magnetic field in our model is presented 
in the magnetic field dependence of the activation energy 
A(-By) = A(n— n c (B»)). We fitted the curves to the data 
using three adjustable parameters po, p\ and A; their nu- 
merical values are indicated on the figure. Whereas there 
is a certain correlation between po and pi, we focus on 
the activation energy A which is well defined and almost 
independent of other parameters. 
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a) Resistivity of Si-MOS sample vs tem- 
zero field and at Bu = IT. Density 



FIG. 9. 

perature at 

n = 0.75 x lO^cm -2 . Reproduced from Fig. 9 of Ref. JsJ. 
Dashed and dash-dotted curves are fitting dependences. 

The carrier density n = 0.75 x 10 11 cm~ 2 is set equal 
to the critical density n = n c (B = 0). By defi- 
nition, the almost horizontal curve for B = corre- 
sponds to A = 0. Fitting the upper curve we ob- 
tained A f it = 0.8 K which is reasonably consistent with 
a value calculated on the basis of our results from Fig. 5: 
A calc = (dA/dn)[n - n c (B = IT)] = 0.6K. For this 
estimate we used the coefficient dA/dn — (4 — 5) [K 



10~ n cm 2 ] measured earlier for variety of high mo- 
bility samples. The lowest temperature point deviates 
substantially from the fitting curve; we presume this is 
due to either electron overheating by ~ 0.07 K [p| or by 
a crossover to another T— dependence. We also didn't 
care about a minor discrepancy for T > 1.5 K; this may 
be eliminated by a simple modification of the model as 
demonstrated further. 

Anyhow, the increase in the measured resistivity by 
four orders of magnitude (which authors of Ref. |34|] 
considered as an extraordinary effect) is at least 10 times 
smaller than the increase demonstrated by the conven- 
tional exponentional dependence Eq. (1) (dashed curve). 
In general, the model curve fits the data rather well and 
some further improvement may be obtained by varying 
the critical index; an example is illustrated in Fig. 8 by 
the dash-dotted curve and the lower formula. The suc- 
cess of this fitting, as well as the consistency between 
the calculated and fitted values for A confirm that the 
action of the magnetic field is simply described by the 
increase in the critical density value n c (B) related to the 
magnetic field induced disorder. 



B. Example 2 

We now turn to the data by Simonian et al. jl5) and 
reproduce them in Fig. 10 a; it shows the temperature de- 
pendence of p for Si-MOS sample measured in fixed var- 
ious in-plane fields between and 1.4 T. The data were 
discussed in Ref. j}4j as a demonstration of an abrupt 
onset of the metallic behaviour and abrupt development 
of the magnetoresistance. 

We note first that the overall p(T) behavior for this 
sample is different from that shown in Fig. |9| here is no 
room for a temperature independent curve in the range 
2 K> T > 0.2 K and the 'separatrix' p(T, n = n c ) is tilted 
[|53[ . Correspondingly, the resistivity curves for this sam- 
ple plotted vs density, for various temperatures (in the 
range (0.2-2) K) would not cross each other at a single 
density |p4ft . 

The carrier density n — 0.883 x 10 n cm~ 2 in Fig. ||a 
was set close to n c so that the effect of the field is 
large. Initially, at B\\ = 0, the resistivity demonstrates a 
'metallic' temperature dependence typical for high mo- 
bility Si-MOS samples (cf. Fig. la), but as field in- 
creases, it transforms to a typical 'insulating' one. We 
estimate the transition between 'metallic' and 'insulat- 
ing' behaviour to occur in the region marked with a 
horizontal arrow, n c (B — 0.93 T) = n. From this 
number and using the empirical dependence n c (B) from 
Fig. 5 we estimated the critical density at zero field, 
n c (B = 0) = 0.76 x 10 11 cm -2 ; this agrees within 5% 
with the number (0.802) given in Ref. @. 

We fitted all curves first with the exponential depen- 
dence p = po + /9iexp(— ), same as Eq. (|^). For the 
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'metallic' region (curves 1, 2, 3), this empirical depen- 
dence with A < is known to describe p(T) for dif- 
ferent samples J3^fi^ i"3| . For the 'insulating' region 
(curves 4, 5, 6), Eq. (|) with A > has a meaning of 
the temperature-activated dependence (though the same 
might be done with a variable-range hopping exponent). 
The three adjustable parameters, A, po and pi, where 
fitted for each curve; their values are given in Table 2. 
We found the fit is rather good for low temperatures, 
T < 2K. 
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FIG. 10. a) Resistivity vs temperature for different 
fixed parallel fields (shown at the top) . Dashed lines show 
fitting curves in the insulating range (curves 1, 2, 3) and 
metallic range (curves 4, 5, 6). Horizontal arrow marks 
an estimated location of the 'critical' trace, b) Two fit- 
ting parameters vs parallel field. Experimental data are 
reproduced from Ref. [[L5| 

Eq. (1) can not obviously describe the high tempera- 
ture region T > 2 K, therefore we also applied a modified 
expression which takes into account the non-exponential 
'tilt' of the experimental data: 



P = [ Po + Pi ex P( 




a- 



T 




(3) 



The two additional parameters a = 0.194 and f3 = 
0.46 — 0.66 are uncorrelated with the other ones, and 
are directly obtained from the tilt at high temperatures. 

Figure [l0| shows that Eq. (||) fits the data rather well 
over the whole range of T whereas for low temperatures, 
T < 2K, both fits, with Eqs. (g) and (g), are almost 
indistinguishable. The resulting parameters for both fits 
are given in Tables 2, and 3, correspondingly. Compar- 
ing the fitting curves with data, we conclude that the 
"abrupt changes" discussed in Ref. are nothing more 
than just a conventional exponential dependence. 

The field dependence of the two principle parameters, 
po and pi is shown in Fig. [l^b: po increases with field 
demonstrating the disordering role of the field (though we 
don't think that the almost linear increase is a universal 
feature for different samples). The strength of the 'resis- 
tance drop' in the 'metallic' range, (po + pi)/ po, decays 
with field, also similar to the decay induced by disorder 
[^,^5|. The critical resistivity, p c , in our model is given 
by the sum po + px; we find a qualitative similarity be- 
tween the two independently calculated curves, p c (B) in 
Fig. 5 and p a + p\ vs B in Fig. [To|b. 



TABLE III. Parameters of the six curves for the fit with 
Eq. (1). n c is in 10 11 cm -2 , po,Pi are in h/e 2 , and A is in 
K. 



curve 


B (T) 


n c (B) 


PO 


Pi 


A flt 


^calc 


1 


1.4 


0.945 


1.45 


0.24 


1 


0.31 


2 


1.2 


0.919 


1.29 


0.35 


0.6 


0.18 


3 


1.0 


0.893 


1.09 


0.6 


0.2 


0.05 


4 


0.9 


0.881 


1.0 


0.8 


-0.15 


-0.01 


5 


0.7 


0.855 


0.9 


1.1 


-0.55 


-0.14 


6 





0.765 


0.45 


2 


-0.71 


-0.59 



TABLE IV. Parameters of the six curves shown in 
Fig. for the fit with Eq. (2). 



curve 


B (T) 


A flt 


Po 


Pi 


1 


1.4 


1.0 


1.82 


0.15 


2 


1.2 


0.6 


1.76 


0.2 


3 


1.0 


0.2 


1.1 


0.6 


4 


0.9 


-0.75 


1.4 


0.95 


5 


0.7 


-0.76 


0.97 


1.47 


6 





-0.75 


0.45 


2.23 



For the curve 3, the fit is not perfect because it oscil- 
lates vs temperature, which is beyond the frameworks of 
the model. The only steep feature may be found in the 
field dependence of p\ at the transition from insulating 
to metallic behavior; this however may be caused by a 
minor mismatching of the oversimplified model curves in 
the metallic and insulating ranges. 
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The 'activation energy' A depends on the choice of the 
model and on average is 2-3 times larger than the values 
estimated from Fig. 5 (cf. the right columns in Table 2). 
This discrepancy is presumably caused by the admixture 
to the data of a non-exponential temperature dependence 
which is essentially strong for the considered sample and 
which can not be separated at low temperatures from the 
exponential one. For example, we were able to achieve a 
reasonable fit of the data with Eq. (2), using the calcu- 
lated in Table 2 values A ca j c = (dA/dB) x (n c — n) when 
we modeled the 'tilted separatrix' in Eq. (^) by a power 
low factor, T~ p , rather than by the polynomial factor. 

In Ref. it was stated that "the effect of magnetic 
field cannot be ascribed solely to a field-induced change 
in the critical electron density" . We verified this possibil- 
ity and arrive at the opposite conclusion. The success of 
our fitting confirms that the action of the parallel field, 
at least, to the first approximation consists in progres- 
sive increase of disorder. This is described by increasing 
the T-independent scattering rate po(B»), decreasing the 
magnitude of the resistivity drop, (pi + po)/po, and in- 
creasing the critical carrier density n c (Bn); the three pa- 
rameters are not independent and within the same mate- 
rial system may be reduced only to the single parameter, 
e.g. n c (B||). 



data |30|,[56| indicates that the mobility in different spin 
and valley sub-bands of conducting electrons are almost 
equal and the intersubband scattering does not play a 
major role in the themperature- and parallel field- de- 
pendence of the resistance. Applying the model of a field 
dependent disorder, we find a qualitative explanation of 
the whole set of presented results. 
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IV. SUMMARY 

To summarize, the results reported in this paper 
demonstrate that the effect of the parallel magnetic field 
on resistivity in high mobility Si-MOS samples (though 
via Zeeman coupling) is similar to that of disorder and, 
to some extent, to that of temperature. In other words, 
parallel magnetic field increases disorder. The temper- 
ature dependence of resistivity for various disorder and 
magnetic fields may be reduced to the dependences of 
the 'critical' density n c on magnetic field and on disor- 
der; the changes in the resistivity drop with temperature, 
(po + pi)/po, caused by disorder and magnetic field may 
be also mapped onto each other. We find a similarity be- 
tween the action of the parallel field and temperature for 
the region of carrier density and field where conduction 
remains 'non-activated'. These findings set constraints 
on the choice of the developed microscopic models. 

The analogy between the parallel field and disorder 
points to the existence of a sub-band of localized car- 
riers. Such band of localized carriers should lift in en- 
ergy as a function of the parallel field, gets 'undressed' 
after passing through the Fermi energy and cause the 
increase in the scattering rate. The strong evidence for 
this mechanism is the variation in the Hall voltage ob- 
served in the vicinity of the metal-insulator transition. 
We can not exclude also the contribution of the inter 
spin-subband scattering, as a complementary feature to 
the above mechanism. However, our Shubnikov-de Haas 
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